Diversity Arrays: a solid state technology for sequence information independent genotyping

Here we present the successful application of the microarray technology platform to the analysis of DNA polymorphisms. Using the rice genome as a model, we demonstrate the potential of a high-throughput genome analysis method called Diversity Array Technology, DArT‘. In the format presented here the technology is assaying for the presence (or amount) of a specific DNA fragment in a representation derived from the total genomic DNA of an organism or a population of organisms. Two different approaches are presented: the first involves contrasting two representations on a single array while the second involves contrasting a representation with a reference DNA fragment common to all elements of the array. The Diversity Panels created using this method allow genetic fingerprinting of any organism or group of organisms belonging to the gene pool from which the panel was developed. Diversity Arrays enable rapid and economical application of a highly parallel, solid-state genotyping technology to any genome or complex genomic mixtures.

PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

DNA sequencing and genotyping by transcriptional synthesis of chain-terminated RNA ladders and MALDI-TOF mass spectrometry

Sets of RNA ladders can be synthesized by transcription of a bacteriophage-encoded RNA polymerase using 3′-deoxynucleotides as chain terminators. These ladders can be used for sequencing of DNA. Using a nicked form of phage SP6 RNA polymerase in this study substantially enhanced yields of transcriptional sequencing ladders. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) of chain-terminated RNA ladders allowed DNA sequence determination of up to 56 nt. It is also demonstrated that A→G and C→T variations in heterozygous and homozygous samples can be unambiguously identified by the mass spectrometric analysis. As a step towards single-tube sequencing reactions, α-thiotriphosphate nucleotide analogs were used to overcome problems caused by chain terminator-independent, premature termination and by the small mass difference between natural pyrimidine nucleotides.

Mouse mammary tumor virus/v-Ha-ras transgene-induced mammary tumors exhibit strain-specific allelic loss on mouse chromosome 4

Hybrid mice carrying oncogenic transgenes afford powerful systems for investigating loss of heterozygosity (LOH) in tumors. Here, we apply this approach to a neoplasm of key importance in human medicine: mammary carcinoma. We performed a whole genome search for LOH using the mouse mammary tumor virus/v-Ha-ras mammary carcinoma model in female (FVB/N × Mus musculus castaneus)F1 mice. Mammary tumors developed as expected, as well as a few tumors of a second type (uterine leiomyosarcoma) not previously associated with this transgene. Genotyping of 94 anatomically independent tumors revealed high-frequency LOH (≈38%) for markers on chromosome 4. A marked allelic bias was observed, with M. musculus castaneus alleles almost exclusively being lost. No evidence of genomic imprinting effects was noted. These data point to the presence of a tumor suppressor gene(s) on mouse chromosome 4 involved in mammary carcinogenesis induced by mutant H-ras expression, and for which a significant functional difference may exist between the M. musculus castaneus and FVB/N alleles. Provisional subchromosomal localization of this gene, designated Loh-3, can be made to a distal segment having syntenic correspondence to human chromosome 1p; LOH in this latter region is observed in several human malignancies, including breast cancers. Evidence was also obtained for a possible second locus associated with LOH with less marked allele bias on proximal chromosome 4.

Typing of urinary JC virus DNA offers a novel means of tracing human migrations

Although polyomavirus JC (JCV) is the proven pathogen of progressive multifocal leukoencephalopathy, the fatal demyelinating disease, this virus is ubiquitous as a usually harmless symbiote among human beings. JCV propagates in the adult kidney and excretes its progeny in urine, from which JCV DNA can readily be recovered. The main mode of transmission of JCV is from parents to children through long cohabitation. In this study, we collected a substantial number of urine samples from native inhabitants of 34 countries in Europe, Africa, and Asia. A 610-bp segment of JCV DNA was amplified from each urine sample, and its DNA sequence was determined. A worldwide phylogenetic tree subsequently constructed revealed the presence of nine subtypes including minor ones. Five subtypes (EU, Af2, B1, SC, and CY) occupied rather large territories that overlapped with each other at their boundaries. The entire Europe, northern Africa, and western Asia were the domain of EU, whereas the domain of Af2 included nearly all of Africa and southwestern Asia all the way to the northeastern edge of India. Partially overlapping domains in Asia were occupied by subtypes B1, SC, and CY. Of particular interest was the recovery of JCV subtypes in a pocket or pockets that were separated by great geographic distances from the main domains of those subtypes. Certain of these pockets can readily be explained by recent migrations of human populations carrying these subtypes. Overall, it appears that JCV genotyping promises to reveal previously unknown human migration routes: ancient as well as recent.

Octamer-based genome scanning distinguishes a unique subpopulation of Escherichia coli O157:H7 strains in cattle

Multilocus-genotyping methods have shown that Escherichia coli O157:H7 is a geographically disseminated clone. However, high-resolution methods such as pulse-field gel electrophoresis demonstrate significant genomic diversity among different isolates. To assess the genetic relationship of human and bovine isolates of E. coli O157:H7 in detail, we have developed an octamer-based genome-scanning methodology, which compares the distance between over-represented, strand-biased octamers that occur in the genome. Comparison of octamer-based genome-scanning products derived from >1 megabase of the genome demonstrated the existence of two distinct lineages of E. coli O157:H7 that are disseminated within the United States. Human and bovine isolates are nonrandomly distributed among the lineages, suggesting that one of these lineages may be less virulent for humans or may not be efficiently transmitted to humans from bovine sources. Restriction fragment length polymorphism analysis with lambdoid phage genomes indicates that phage-mediated events are associated with divergence of the lineages, thereby providing one explanation for the degree of diversity that is observed among E. coli O157:H7 by other molecular-fingerprinting methods.

Asian genotypes of JC virus in Native Americans and in a Pacific Island population: Markers of viral evolution and human migration

The human polyomavirus JC (JCV) causes the central nervous system demyelinating disease progressive multifocal leukoencephalopathy. Previously, we showed that 40% of Caucasians in the United States excrete JCV in the urine as detected by PCR. We have now studied 68 Navaho from New Mexico, 25 Flathead from Montana, and 29 Chamorro from Guam. By using PCR amplification of a fragment of the VP1 gene, JCV DNA was detected in the urine of 45 (66%) Navaho, 14 (56%) Flathead, and 20 (69%) Chamorro. Genotyping of viral DNAs in these cohorts by cycle sequencing showed predominantly type 2 (Asian), rather than type 1 (European). Type 1 is the major type in the United States and Hungary. Type 2 can be further subdivided into 2A, 2B, and 2C. Type 2A is found in China and Japan. Type 2B is a subtype related to the East Asian type, and is now found in Europe and the United States. The large majority (56–89%) of strains excreted by Native Americans and Pacific Islanders were the type 2A subtype, consistent with the origin of these strains in Asia. These findings indicate that JCV infection of Native Americans predates contact with Europeans, and likely predates migration of Amerind ancestors across the Bering land bridge around 12,000–30,000 years ago. If JCV had already differentiated into stable modern genotypes and subtypes prior to first settlement, the origin of JCV in humans may date from 50,000 to 100,000 years ago or more. We conclude that JCV may have coevolved with the human species, and that it provides a convenient marker for human migrations in both prehistoric and modern times.

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6–40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

A novel modifier gene for plasma von Willebrand factor level maps to distal mouse chromosome 11

Type 1 von Willebrand disease (VWD), characterized by reduced levels of plasma von Willebrand factor (VWF), is the most common inherited bleeding disorder in humans. Penetrance of VWD is incomplete, and expression of the bleeding phenotype is highly variable. In addition, plasma VWF levels vary widely among normal individuals. To identify genes that influence VWF level, we analyzed a genetic cross between RIIIS/J and CASA/Rk, two strains of mice that exhibit a 20-fold difference in plasma VWF level. DNA samples from F2 progeny demonstrating either extremely high or extremely low plasma VWF levels were pooled and genotyped for 41 markers spanning the autosomal genome. A novel locus accounting for 63% of the total variance in VWF level was mapped to distal mouse chromosome 11, which is distinct from the murine Vwf locus on chromosome 6. We designated this locus Mvwf for “modifier of VWF.” Additional genotyping of as many as 2407 meioses established a high resolution genetic map with gene order Cola1-Itg3a-Ngfr-Mvwf/Gip-Hoxb9-Hoxb1-Cbx·rs2-Cox5a-Gfap. The Mvwf candidate interval between Ngfr and Hoxb9 is ≈0.5 centimorgan (cM). These results demonstrate that a single dominant gene accounts for the low VWF phenotype of RIIIS/J mice in crosses with several other strains. The pattern of inheritance suggests a gain-of-function mutation in a unique component of VWF biosynthesis or processing. Characterization of the human homologue for Mvwf may have relevance for a subset of type 1 VWD cases and may define an important genetic factor modifying penetrance and expression of mutations at the VWF locus.

Direct genetic analysis by matrix-assisted laser desorption/ionization mass spectrometry

An approach to analyzing single-nucleotide polymorphisms (SNPs) found in the human genome has been developed that couples a recently developed invasive cleavage assay for nucleic acids with detection by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The invasive cleavage assay is a signal amplification method that enables the analysis of SNPs by MALDI-TOF MS directly from human genomic DNA without the need for initial target amplification by PCR. The results presented here show the successful genotyping by this approach of twelve SNPs located randomly throughout the human genome. Conventional Sanger sequencing of these SNP positions confirmed the accuracy of the MALDI-TOF MS analysis results. The ability to unambiguously detect both homozygous and heterozygous genotypes is clearly demonstrated. The elimination of the need for target amplification by PCR, combined with the inherently rapid and accurate nature of detection by MALDI-TOF MS, gives this approach unique and significant advantages in the high-throughput genotyping of large numbers of SNPs, useful for locating, identifying, and characterizing the function of specific genes.

Database resources of the National Center for Biotechnology Information

In addition to maintaining the GenBank® nucleic acid sequence database, the National Center for Biotechnology Information (NCBI) provides data analysis and retrieval resources that operate on the data in GenBank and a variety of other biological data made available through NCBI’s Web site. NCBI data retrieval resources include Entrez, PubMed, LocusLink and the Taxonomy Browser. Data analysis resources include BLAST, Electronic PCR, OrfFinder, RefSeq, UniGene, HomoloGene, Database of Single Nucleotide Polymorphisms (dbSNP), Human Genome Sequencing, Human MapViewer, GeneMap’99, Human–Mouse Homology Map, Cancer Chromosome Aberration Project (CCAP), Entrez Genomes, Clusters of Orthologous Groups (COGs) database, Retroviral Genotyping Tools, Cancer Genome Anatomy Project (CGAP), SAGEmap, Gene Expression Omnibus (GEO), Online Mendelian Inheri­tance in Man (OMIM), the Molecular Modeling Database (MMDB) and the Conserved Domain Database (CDD). Augmenting many of the Web applications are custom implementations of the BLAST program optimized to search specialized data sets. All of the resources can be accessed through the NCBI home page at: http://www.ncbi.nlm.nih.gov.

A single nucleotide polymorphism in the RAD51 gene modifies cancer risk in BRCA2 but not BRCA1 carriers

BRCA1 and BRCA2 carriers are at increased risk for both breast and ovarian cancer, but estimates of lifetime risk vary widely, suggesting their penetrance is modified by other genetic and/or environmental factors. The BRCA1 and BRCA2 proteins function in DNA repair in conjunction with RAD51. A preliminary report suggested that a single nucleotide polymorphism in the 5′ untranslated region of RAD51 (135C/G) increases breast cancer risk in BRCA1 and BRCA2 carriers. To investigate this effect we studied 257 female Ashkenazi Jewish carriers of one of the common BRCA1 (185delAG, 5382insC) or BRCA2 (6174delT) mutations. Of this group, 164 were affected with breast and/or ovarian cancer and 93 were unaffected. RAD51 genotyping was performed on all subjects. Among BRCA1 carriers, RAD51-135C frequency was similar in healthy and affected women [6.1% (3 of 49) and 9.9% (12 of 121), respectively], and RAD-135C did not influence age of cancer diagnosis [Hazard ratio (HR) = 1.18 for disease in RAD51-135C heterozygotes, not significant]. However, in BRCA2 carriers, RAD51-135C heterozygote frequency in affected women was 17.4% (8 of 46) compared with 4.9% (2 of 41) in unaffected women (P = 0.07). Survival analysis in BRCA2 carriers showed RAD51-135C increased risk of breast and/or ovarian cancer with an HR of 4.0 [95% confidence interval 1.6–9.8, P = 0.003]. This effect was largely due to increased breast cancer risk with an HR of 3.46 (95% confidence interval 1.3–9.2, P = 0.01) for breast cancer in BRCA2 carriers who were RAD51-135C heterozygotes. RAD51 status did not affect ovarian cancer risk. These results show RAD51-135C is a clinically significant modifier of BRCA2 penetrance, specifically in raising breast cancer risk at younger ages.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.